Localization of gravity in brane world cosmologies 
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The most remarkable and interesting feature of brane world scenario is the use of bulk's curvature 
to localize gravity on the brane, albeit with fine tuning of the brane and bulk parameters. For FRW 
expanding universe on the brane, it is a moving hypersurface in Schwarzschild anti de Sitter bulk 
spacetime. We show that zero mass gravitons have bound state on the brane for suitable values of 
brane and bulk parameters. There occur various cases giving rise to different cosmological models, 
in particular we discuss a model with positive cosmological constant on the brane. 
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The view that our Universe might actually have dimen- 
sions more than four is anchored on the recent develop- 
ments in string and M-theories in which gravity arises as 
a truly higher dimensional interaction. Only in the low 
energy limit it manifests in the familiar 4-dimensional 
general relativity (GR). This has initiated a lot of activ- 
ity in recent times. Though the basic idea was already 
there in the form of Kaluza-Klein (KK) theories, the re- 
cent spurt is primarily due to the possibility which helps 
solve the mass hierarchy problem in the standard model 
of particle physics. 

In some of these models [1] the extra dimensions can 
be as large as millimeter which is however less than the 
current observational limits on low scale gravity. No- 
table are the models in which the extra dimensions are 
allowed to be of infinite extent [2, 3]. These models have 
Z2-symmetry, which is motivated by the reduction of M 
theory to x Eg, heterotic string theory [4] . The single 
brane Randall-Sundrum (RS) model [3] has attracted a 
lot of interest and activity. In this model the Minkowski 
fiat brane in 5-D anti de Sitter (AdS) bulk has a positive 
tension. It is then possible to recover Newton's inverse 
square law with correction term which arises from 
massive KK modes contribution. There have been vari- 
ous generalizations of this [5] in the form of thick branes 
[6] , AdS branes [7] and brane models without Z2 symme- 
try [8]. 

In the overall view of the brane world scenario, all mat- 
ter fields live on the 3-brane which is the 4-D physical 
Universe while gravity can propagate in the extra dimen- 
sions, say 5-D bulk. The bulk and brane spacetimes are 
joined together through the Israel junction conditions [9]. 
Since, the standard Einstein equations on the brane are 
modified by the bulk effects it opens a whole new vista 
for investigation of astrophysical and cosmological con- 
sequences of these models (see for eg. [10]). The con- 
nection with CFT correspondence has also been studied 
(see for eg. [11]). For complete solution of the problem, 
one has to solve the A-vacuum equation in 5-D for the 



*e-mail: param@iucaa.ernet.in 
te-mail: nkd@iucaa.ernet.in 



bulk spacetime and simultaneously the modified Einstein 
equation which in addition to the square of stress tensor 
also contains the projection of the bulk Weyl curvature 
on the brane (dark radiation) and then the two solutions 
should be joined together with the Israel junction con- 
ditions [9]. It is by all means a very formidable task 
and it is therefore no surprise that there exist only few 
complete solutions to date. The two important cases are 
fiat/ vacuum brane with AdS bulk and FRW brane with 
Schwarzschild- AdS (S-AdS) bulk [12-14]. 

The S-AdS bulk-brane system is composed of two 
patches of S-AdS bulk having in general different mass 
parameters for the black hole with the brane located at 
y = yir), where y(r) is determined by solving the Is- 
rael junction conditions [13]. From the resulting equa- 
tion for brane trajectory one finds that in S-AdS bulk, 
the FRW brane will in general be moving unless the pa- 
rameters are properly fine tuned (eg. for the RS case 
cr = (3/47rG5Z), A4 — 0, where a is the brane tension, 
G5 is the 5-D gravitational constant and I is the radius 
of curvature of the bulk spacetime). The extra dimen- 
sion is a radial coordinate of the bulk and imposing Z2 
symmetry across the brane demands that the mass pa- 
rameters of both the patches to be same. By fine tuning 
parameters one can obtain static branes which cannot 
harbour expansion which is essential for realistic cosmo- 
logical models. A slight off tuned value of a or non zero 
value of black hole mass would set it moving. 

Though localization of gravity for the AdS bulk with 
flat brane and for some generalizations of it has been 
established [7, 15, 16], it has not yet been investigated 
for the S-AdS bulk with FRW brane. The problem be- 
comes difficult if one notes that unlike the RS case, the 
brane is dynamic in the static bulk and it is non triv- 
ial to fix the boundary conditions on the modes. For 
localizability of gravity on the brane, first there should 
exist bound normalizable mode for zero mass graviton 
and plus there should also exist KK modes to give the 
correction to the Newtonian gravity. All this could be 
studied by considering perturbation of the bulk metric 
and the brane motion determined by the Israel junction 
conditions. This is what we shall concern ourselves in this 
investigation which should be quite indicative of localiza- 
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tion of gravity on the branc. For actual demonstration 
of recovering Newtonian gravity with the correction, one 
has to do propagator analysis which is hard to carry out 
for a non-static brane and we shall not attempt that here. 

It is well known that localizability is very sensitive to 
fine tuning of parameters. A priori, there is no well es- 
tablished criterion to check this. For instance there does 
exist a bulk spacetime, which is an exact solution of the 
A- vacuum equation, for which gravity is non localizable 
on the brane [17]. This is the case of Nariai metric which 
has non zero Weyl curvature. Note that Weyl is non zero 
for S-AdS as well. It is therefore pertinent to find un- 
der what conditions do zero mass gravitons have bound 
state on FRW brane? This is the most critical question 
for brane world cosmologies and that is what we wish to 
address in this letter. 

Since our brane is a hypersurface around the black 
hole, any movement of the brane towards or away from 
the black hole would be interpreted as contraction or ex- 
pansion by the observers on the brane. That is how ex- 
pansion is directly related to motion of the brane. To 
see whether zero mass graviton has a bound normaliz- 
able state in these cosmological models one has to solve 
for the graviton perturbation equation by taking into ac- 
count the location of the brane. Unlike the static cases 
(like the RS case) the position of the FRW brane would 
not be fixed in the bulk and one would have to take 
into account the trajectory of the brane to understand 
the localizability. The brane trajectory found through 
Israel matching conditions would be given by a Fried- 
mann like equation which determines the position of the 
brane. In the following we shall hence study the metric 
perturbations of the S-AdS bulk and obtain a potential 
in the graviton wave equation which would determine 
the fate of localizability once the brane trajectory (and 
hence the location of the brane) is taken care of. We 
shall then show how the RS case can be recovered in this 
scenario, then we shall in particular consider the case of 
fc = 1,A4 > 0, while a comprehensive discussion of all 
possible cases would be done elsewhere [18]. 

In the five dimensional bulk we have the S-AdS metric, 
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Here k = 0, ±1 is the curvature index and M is the mass 
parameter of the bulk black hole. 

The static limit would be given by goo = which leads 
to. 
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and since the metric is spherically symmetric this also de- 
fines the location of event horizon of the spacetime. The 



above metric is the solution of the A-vacuum equation, 

Gab = -A5gab, A5 = -6/l^. (4) 

The Latin indices which label the bulk spacetime (.T'^,y) 
run from 0...4 and the Greek indices labelling brane 
spacetime (x^) would run from 0...3. We consider the 
metric perturbations of the above metric , i.e. hab = 

9ab — Sab ■ would take the metric perturbations in 
the extra dimension to vanish, hty = hry = hgy = 
hffiy = hyy = 0. Wb wouM furthcr impose the transverse- 



traceless gauge conditions, V^h/j,,, = 0, /i''^ 
hence the wave equation turns out to be. 



and 



Ri^^ha' = A,hab. (5) 



Assuming that wavelength of the gravitons is much 
smaller than the radius of curvature of the background 
and their amplitude is very small we can work under lin- 
earized approximation. We can further choose a constant 
vector field which satisfies habU^ = and hence we are 
effectively left with two independent modes. We solve 
this wave equation with the ansatz that habix^,y) = 
hab{x'^) ^{y)- Substituting this in eq.(5) and using to^ as 
the constant of separation of variables, the y dependence 
turns out to be 
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where prime denotes a derivative with respect to y. This 
equation can be written down in the form 
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which can be transformed into a wave equation form with 
the transformation *(j/) = 4'{y)ij{y) where 



<p{y) = cexp(^-^ J ai{y) dy^ . 



(8) 



Here c is a constant and yi and yf belong to the interval 
over which 01,02 are continuous and Oi has a continu- 
ous derivative. This transformation eliminates the first 
order derivative term in eq.(6) and we get the desired 
Schroedinger like equation 
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The potential consists of two parts, V = Vf + Vm, which 
are given by 
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The presence of interaction part is peculiar to the coor- 
dinate system we are working with. Obviously, the shape 
of the potential will play determining role in localization 
of gravitons. The potential V for all the cases provides 
a negative infinite well at the event horizon of the black 
hole of mass M, which holds irrespective of the matter 
distribution on the brane. The potential for the massless 
and massive modes for the case A; = 1, M 7^ is shown 
in Fig. 1. The potential for the massless mode can be 
viewed as the m — * limit of the potential for massive 
modes for which ^ > for large y. Such a behaviour is 
a generic feature of various cosmological scenarios. 

The Israel junction conditions determine the relation- 
ship between the bulk and brane parameters, 

A4 - — + —^G^a , G4 = —G^cr (12) 

and also the motion of the brane given by the equation 
[13], 

Here dot refers to derivative relative to proper time t 
and p is the energy density on the brane and the critical 
tension cro = 3/47rG5/, which determines the sign of A4. 
The y coordinate now plays a dual role. It not only tells 
us about the effect on graviton fiuctuations due to extra 
dimension but also parameterizes the brane trajectory. 
The metric on the brane is FRW and is given by, 

dsl = -dr'+y\r)[^^ + dnl^ (14) 

where dr^ = ea:p(2/3(t))(l - exp{-Ap{t)){dy / dtf)dt^ . 
The physically interesting branes would lie outside the 
horizon which would also serve as the high energy cut- 
off. For localization of massless mode we require it to be 
bounded as well as normalizable. The form of the poten- 
tial for various choice of cosmological parameters is such 
that boundedness of massless mode is guaranteed once 
the brane crosses horizon. However, the normalizability 
of the massless mode depends on the asymptotic behav- 
ior of the potential and the form of the potential reflects 
that branes which are either static or ever expanding in 
the bulk would harbour localization. 

As is clear from eq. (13) that motion of brane is deter- 
mined by both black hole mass and energy distribution 
on the brane. The localization would therefore depend 
on energy distribution in both bulk and brane. The lo- 
calization of the massless mode depends on whether the 
brane expands forever or has a bounce, and it can be 
easily checked from eq.(13) that such a behaviour of the 
brane trajectory is unchanged even if the matter density 
terms are excluded. For simplicity we would henceforth 
exclude matter density terms from our analysis of brane 
trajectory and hence the brane trajectory equation be- 
comes 

y'-^y'- — + fc = o. (15) 



Note that this equation is non linear in the highest order 
of derivative and hence it will not have unique solution. 

It should be noted for dynamic branes the proper time 
as measured by the observers on the brane is not the 
same as the time component of the bulk metric. For 
the observer on the brane to interpret the massless mode 
correctly one thus has to assume that y(t) is a slowly 
varying function. Further, fixing boundary conditions 
for the modes on a moving boundary is a very difficult 
task. However, in a cosmological setting of FRW branes 
we can safely assume that brane is slowly expanding and 
these problems can be bypassed. Our results for dynamic 
branes would hence hold in the above cosmological sce- 
nario. 
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FIG. 1: Potential plot for M / and = 1. The dashed and 
dark lines indicate potential for massless and massive mode 
respectively. 

Now we would see how the above metric perturba- 
tion equation and the brane dynamics equation yield RS 
model. For this fc = M = 0, cr = cro, the metric eq.(l) 
takes the form 

ds^ = -^dt^ + !^dy^+y^ {dr^ + r^dr^ + sin^ Od^^) . 

' y 

(16) 

Using the transformation r] = I \n{l/y) one can rewrite 
this metric in the familiar RS form. Note that j/ is a radial 
coordinate of the black hole in the bulk and it ranges 
from to 00. The above transformation would yield r] 
ranging from —00 to 00. Imposition of Z2 symmetry 
across the brane (i.e. matching the extrinsic curvatures 
on both sides of the brane) does not give any condition 
on y coordinate, however in r] coordinate it demands ??(_) 
to be identified with ?7(_|_) and hence the metric contains 
1 77 1 in the warp factor, 

ds^ = dri^ + e-2|''l/'(-dt2 + dxl + dxl + dxl). (17) 

RS brane is static and located at 77 = 0, which means 
Ski y = I. Note that the shape and form of the poten- 
tial in the Schroedinger equation critically depends on 
the coordinates one is working with. In the 77 coordi- 
nates it is the presence of jr;! which gives a Dirac delta 
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potential at the location of the brane. However, work- 
ing with the y coordinate one gets from eqs.(lO) & (H), 
V = -1/Sy'^ + m^(l - P/2y'^). At the location of the 
brane the massless mode is bound and the asymptotic 
behaviour of the potential suggests that it is normaliz- 
able. The massive modes for which m > 1/21, > 0, 
are unbound. Solving near the brane one can obtain the 
wavefunctions from eq.(9) which turn out to be similar to 
those obtained by RS, with the required ml suppression 
of KK modes on the brane [18]. The fact that V turns 
positive for massive modes gives the continuum spectrum 
and the RS correction to Newtonian potential. The form 
of the potential also suggests that there exist discrete 
modes for m < 1/21. The appearance of discrete modes 
in this case is due to our choice of the coordinates. This 
is clear from the fact that unlike RS coordinates our ex- 
tra dimension is a radial coordinate in the bulk and our 
brane is curved. However, it should be noted that the cor- 
rection to Newtonian gravity is not expected to change 
since the Green's hmction used by Garriga and Tanaka 
[19] to evaluate the force law between two point sources 
on the brane in the RS case is inert under our transfor- 
mation of coordinates from y to rj. In the y coordinates 
their Green's function can be written as 



Gr{x,x') = 



Um{y)um{y') dm 



[w + ie)^) 
(18) 



where Umiv) are the wavefunctions for massive modes. 

In the stationary case the Green's function between two 
points at the location of the brane, i.e. y — y' = I yields, 



G{x\l,x'\l) 



1 



(19) 



which loads to the same I'^/r^ correction to Newtonian 
gravity. 

For negative A4 brancs the brane trajectory exhibits 
a bounce, like for the case k = — 1, M = 0, A4 < it is 
given by ?/(t) = sin(-y/— A4/3T)/-y/— A4/3. The potential 
is similar as in RS case, apart from that it blows up at 
origin and there is an infinite negative well eX y = I, but 
the brane trajectory is such that it does not yield the re- 
quired behaviour of the potential for the ground state 
wavefunction to be normalizable. Hence the massless 
mode would not be localized for the negative A4 brane 
(but that would be for the positive A4 brane because the 
solution of the brane motion equation does permit nor- 
malized wavefunction). We thus recover the well-known 
result for AdS bulk and negative/positive A4 branes [7]. 
The behavior of brane trajectory also plays critical role 
in localization of gravity on the brane. 

Now we turn to FRW brane. Eq.(15) yields host of 
solutions for interesting cosmological scenarios including 
the inflationary solutions for all values of k with a pos- 
itive A4 on the brane, which is also favored by current 



observations of type la supernovae [20]. A detailed dis- 
cussion of all these cases would be done elsewhere [18], 
here we would as a representative consider the case of 
fc = 1, M ^ and A4 > 0. 

Solving eq.(15) for this case we get, 



1 -|- n sinh(2a;) — cosh(2a;) 



1/2 



(20) 



where n = 2^fMK^/2) and x = ^ ^-a/Zt. The behavior 
of potential V for this case is shown in Fig.l. Note that V 
blows up at T/ = and there is an infinite well at y = y^, 
we shall therefore restrict to y > yh- Since the form of 
the potential for localization requires an ever expanding 
brane, only n > 1 is possible. The brane emerges out of 
the event horizon, expanding from y = at r = like a 
white hole [21] and would expand for ever, exponentially 
for large r. The Hubble parameter for A; = 1, M 7^ 
cosmological model is given by 



H 



n cosh(2a;) — sinh(2a;) 



1 + n sinh(2a;) — cosh(23;) 



(21) 



implying an inflationary universe at large t. 

For the massless mode the potential profile favors 
boundedness and normalizability. The massive modes 
woiild be unbounded and would contribute to correction 
over Newtonian potential. In order to show the correction 
suggested by massive modes we would first parameterize 
the location of the brane off the horizon by a = 1 — M/y"^. 
Solving the Schroedinger equation, eq.(9), in the approx- 
imation M <^ P and near the horizon {yh ^ VM) we 
get, 

V(y) = Vy{Cii-y/2{'^y) + C2i^/2{j^y)) (22) 

where 7 = ^/l^^W{^^^4a){M^^2aP)/W and u = 
■\/(l — a){2aP — M)ml / M . The boundary conditions 
on the wavefunction at the brane in a cosmological setting 
leads to Ci = (mO^^+^^/^MT/* with 7 < 2 and C2 = 0. 
Close to the horizon v ~ ml/\[M and since y ^ \fM, 
the argument in the Bessel function as in the RS case 
is proportional to ml at the location of the brane. This 
suggests a P/r^ correction to Newtonian potential due 
to massive modes. 

The modifications to the standard GR would be most 
prevelant near the event horizon which is the high energy 
end and marks a cut off for the scale factor y(r). As the 
brane moves out and expands, the potential as shown in 
Fig.l becomes shallower and the high energy modifica- 
tions die out with time. In particular, the RS correction 
to the Newtonian potential will die out as the universe 
expands. 

The overall picture that emerges is that we could have 

different FRW cosmological models on the brane which 
are anchored onto S-AdS bulk. The curvature of the 



5 



bulk spacctimc and the branc trajectory are key players 
in the localization of gravity on the branc. The ultimate 
evolutionary fate of models is decided by the black hole 
mass and the cosmological constant on the brane. We 
have thus shown that S-AdS bulk does allow localization 
of gravity on FRW brane with A4 > 0. 

In conclusion, we could say that brane cosmologies are 
firmly anchored in respect of localization of gravity. We 
have shown that the potential in the Schroedinger equa- 
tion and the brane dynamics lead to bound state for zero 



mass graviton on a slowly moving brane in a cosmological 

setting. 
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